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Abstract 

We investigate intercellnlar insnlin dynamics in an array of diffnsively conpled 
pancreatic islet /3-cells. The cells are connected via gap jnnction conpling, where 
nearest neighbor interactions are inclnded. Throngh the mnltiple scale expansion 
in the semi-discrete approximation, we show that the insnlin dynamics can be 
governed by the complex Ginzburg-Landau eqnation. The localized solntions of 
this eqnation are reported. The results suggest from the biophysical point of view 
that the insulin propagates in pancreatic islet /3-cells using both temporal and 
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spatial dimensions in the form of localized modulated waves. 

Keywords: insulin, islet /3-cells, complex Ginzburg-Landau equation, modulated 
waves. 


1. Introduction 


Blood glucose levels are controlled by a complex interaction of multiple chemi¬ 
cals and hormones in the body. The metabolism of glucose into the /3-cells leads 
to increase in the adenosine triphosphate concentration, closure of ATP-sensitive 
channels, depolarization of the /3-cell membrane and opening of the voltage- 
dependant (7^+ channels, thereby allowing influx [l|]. The resultant rise in 
intracellular concentration in the ^d-cell triggers insulin secretion. Insulin, 
which is secreted from pancreatic ^d-cells is the key hormone regulating glucose 
levels. 


Physiological responses generated within a cell can propagate to neighboring cells 
through intercellular communication involving the passage of a molecular signal to 


a bordering cell through a gap junction [2|,l3[, through extracellular communication 
involving the secretion of molecular signals j2| such as hormones, neurotransmit¬ 
ters, etc., and also trough extracellular calcium signaling jj, 5|. Oscillations of 
rather than metabolism in the ^d-cell are thought to be the direct cause of 


these oscillations in insulin secretion j^. Sneyd et al. in proposed a dynamical 
model of such oscillations, which assume gap junction diffusion of inositol-1,4,5- 
triphosphate between adjacent cells. The diffusion of inositol-1,4,5-triphosphate 
between cells then initiate not only oscillations but also insulin oscillations in 
adjacent cells. 


0 


It is well known that the dynamics of the insulin is very relevant because it is 
related to the onset of the pathologies such as diabetes caused by elevated blood 
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glucose. Careful diabetes mellitus self-management is essential in avoiding chronic 


complications that compromise health, and is characterized by many and often not 
readily observable clinical effects js]. Along the same line, attention has been paid 
with the increased emphasis on derangements of the sensitivity of tissues to insulin 
in diverse pathological conditions like diabetes, obesity and cardiovascular diseases 
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ll| . Therefore, there is an urgent need for improved diagnostic methods 
that provide more precise clinical assessments and sensitive detection of symptoms 
at earlier stage of the disease. This may be facilitated by improved mathemat¬ 
ical models and tools related to interrelationship dynamics among physiological 
variables implicated in the glucose-insulin system. This assumption motivates the 
present work, where a diffusive model of coupled pancreatic islet /3-cells is inves¬ 
tigated. The cells are connected via the gap junction coupling which consists of 
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a mechanism used by cells to coordinate and synchronize their information 
We then investigate a clear analytical solution describing the dynamics of insulin 
in the system. By means of the multiple scale expansion in the semi-discrete ap¬ 
proximation, we use the Lienard form of a diffusive system to obtain the complex 
Ginzburg-Landau (CGL) equation, which describes the evolution of modulated 
waves in this system. We obtain an expression of the hormonal wave by using of 


the envelope soliton solution of the GGL by Pereira and Stenflo [1^, and Nozaki 
and Bekki . The solution reveal that the hormonal wave is a localized nonlinear 
solution which propagates in the form of a breather-like coherent structure. 


The rest of the paper is organized as follows: in Section [2], we present the mathe¬ 
matical model. In Section [3l we hnd an envelope soliton in the model by applying 
the multiple scale expansion in the semi-discrete approximation. Our work is 
summarized in Section 01 
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2. Mathematical model 


We propose in this paper a diffusive model of coupled pancreatic islet /9-cells, in 
which the cells are connected via the gap junction coupling. Diffusive cell models 
with gap junction coupling are quite interesting in describing and characterizing 
quasi-perfect intercellular communication [3]. For the mathematical modeling, let 
us consider and yn as the intracellular concentrations of glucose and insulin in 
the nth cell, respectively. The equations of N coupled cells in the system through 
the gap junction are given by 


Xn CLlXfi 0‘2^nyn T 


( 1 ) 


ijn = D{yn +1 - 2yn + yn-l) + hXn - 


( 2 ) 


where n = 1,...,A^. The parameter ai is the rate constant which represents 
insulin-independent glucose disappearance, 02 is the rate constant which repre¬ 
sents insulin-dependent glucose disappearance. In other terms, 02 describes the 
modulation of the effective kinetic constant of the glucose utilization by insulin 
action. The parameter 03 is the glucose infusion rate, bi is the rate constant 
which represents insulin production due to glucose stimulation and 62 is the rate 
constant which represents insulin degradation. The parameter D is the coupling 
strength of the gap junction. We have considered two nearest neighbors coupling 
in a weak coupling regime. A weak coupling between neighboring cells is a situa¬ 
tion that arises in the study of bursting activity in the /3-cell islets of the pancreas, 
which secrete insulin in response to glucose level [l^, 151. The model also predicts 
that oscillations occur if there is sufficient diffusion {D > 0.1) to create adequate 
concentrations mixing in the reacting layers of the cells 16|. 


Note that the model is nonlinear, due to the presence of the bilinear term between 
the internal variable y{t) representing the insulin action and the variable x{t) 
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representing in the first equation the glucose concentration. The parameter values 
chosen are those from the literature of healthy volunteer undergoing the intra 
venous gluc ose tolerance test. In a clinical experiment conducted and reported 


in Ref. 


17|, ten healthy volunteers (5 males and 5 females) participated in the 


study. As indicated by Gaetano and Arino [17[, all subjects had negative family 
and personal histories for diabetes mellitus and other endocrine diseases, were on 
no medications, had no current illness and had maintained a constant body weight 
for the six months preceding the study. In the present paper, we have taken the 
data of the two hrst subjects as listed in the Table 1. 


Subjects 

Ol 

02 

03 

bi 

&2 

1 

0.0226 

3.8 X 10-® 

1.56 

0.0022 

0.0437 

2 

0.0509 

1.29 X 10-^ 

4.02 

0.0096 

0.2062 

Table 1: Parameter values 1 

Zl- 


Many papers have been published detailing the occurrence mode oscillations of 
pancreatic islet /3-cells. In the present paper, the system being nonlinear, we are 
mainly interested by nonlinear solutions that can describe the nonlinear dynamics 
of insulin. It is convenient to transform the system into the wave form. This is 
achieved by differentiating the second equation and substituting into the ob¬ 
tained second-order ordinary differential equation. The above transformations did 
not fundamentally affect the structure of the system, but allow us to conveniently 
write Eqs. ([1]) and ([2]) in a Lienard form that is a second-order differential equation 
with a small damping term. The governing equation for the insulin concentration 
then reads 


ijn + ^‘oVn + + ^iyn)yn + + A' — Do(yn-|-l ~ -|- ^n-l) 

+ Diijjn+l — 2yn + ijn-l) + D2{yn+1 ~ 

where = 01 ^ 2 , = «i + h, I'l = 02 , Ai = 0262 , Dq = aiD, Di = D, D 2 


( 3 ) 


CI2D 
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and K = —bia^. 


For such equations, perturbation approaches are used to obtain nearly exact solu¬ 
tions. Accordingly, we introduce a new variable ipn such that 

yn = e'ljjn- (4) 

As we are interested for solution in a weakly dissipative medium, we assume the 
parameters uq and Di perturbed at the order e^. Keeping the hrst nonlinear term 
of the development, Eq. ([3]) reads 

'tpn + + e^l'0n)'0n + ~ ‘^4’n + '0n-l) 

( 5 ) 

+ Diitjjn+l ~ ‘^'<Pn + t/'n-l) + cF*2('0n+l ~ + V’n-l)'0n- 


The solutions of this equation are regarded up of carrier waves modulated by 
envelope signal, called envelope solitons which appears naturally for most weakly 


dispersive and nonlinear systems in the small amplitude limit 


18|. In the next 


section, the multiple scale expansion in the semi-discrete approximation is used 
to hnd the envelope soliton solution of Eg. (|3]). The method has been found a 


powerful tool in solving similar equation 181. 


3. Multiple scale expansion in the semi-discrete approximation 


Referring to the multiple scale expansion 


19 


20 | , we proceed by making a change 


of variables according to the space and new time scales = e"'z and T„ = 
respectively. That is we are looking for solution y{z, t) depending on these new 
set of variables as a perturbation series of functions 


{Zq, Zi, Z 2 , ..., To, Ti, T 2 , ...), 


( 6 ) 


n=l 


where Zn and T„ are treated as independent variables. 
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3.1. Equation of motion of the amplitude 

The semi-discrete approximation is a perturbation technique in which the carrier 
waves are kept discrete while the amplitude is treated in the continuum limit. For 
this, we look for modulated wave solution of the form 


fjn -f + e[Fo„ + -f + 0{e‘^), (7) 


where 6^ = qn — cot, q is the wave vector and u is the frequency. The procedure 
will consist of replacing the form of the solution Eq. ([7]) and its derivatives in the 
different terms of Eq. (|5]). We then group terms in the same power of e, which 
leads us to a system of equations. Each of those equations will correspond to each 
approximation for specihc harmonics. 


Substituting Eq. ([7]) in Eq. (|5]) gives 

(Fi,n ~ -f e(Fo,n) + — 4za;F2 „ — 

+ + ^Fo,n + ^F2^„e^*®"] -|- e^z/o(Fi,n ~ 

+ eFi,„Fo,„F®'^ + eFo,„(Fi,„ - *a;Fi,„)F®~] + + 2Fi,„F*J 

= +2e2Ai(Fi,„Fo,„ + F*„F2,„)F^"] + Fo(Fi,„+iF'^“ + Fi,„_ie-*‘'“ - 2Fi,„)F^" 

+ eDo{Fo,n+i + Fo,„_i - 2Fo,„) + eFo(F2,„+ie2*''“ + F2,„_ie-2*''“ - 2F2,^)e^^^- 
+ e2Fi[Fi,„+iF^“ -1- Fi,._ie-*'^“ - 2A,n - zn;(Fi,„+iF^“ + Fi,„_ie-*'^“ - 2Fi,„)]F®'^ 

+ 6F2Fi,„(Fi,„+iF''“ + Fi,„_ie-*'^“ - 2Fi,„)e2*^" + e‘^D2Fi^n{Fo,n+i + Fo,n-i - 2Fo,„)F®'^ 
+ e2F2Fo,„(Fi,„+iF''“ + Fi,„_ie-*'^“ - 2Fi,„)F^". 

( 8 ) 


Since the envelope function varies slowly in space and time, we use the continuum 
approximation for F in a multiple scale expansion such that 


Fn±i = F ± e 


dF 


±e^ 


dF 


+ 


d‘^F 


dZ2 2 dZf 


+ 0{e^), 


(9) 
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and 



( 10 ) 


Equating dc, first-, and second-harmonic terms, we get respectively 



( 11 ) 


Ai — iuPi — 4D2 sin^(|) 2 

— Tw:;— ■ A — -E, 


31]2 + l6Dosin^(|) 


( 12 ) 



(13) 


In the above calculation, we have used the dispersion relation for the carrier 
wave 


a;2 = _^4Dosin^(|), 


(14) 


obtained by linearizing Eq.(8). As we observe in Fig. (1), the corresponding linear 
spectrum for the hrst two subjects of Ref. jl^ is related to the system parame¬ 
ters. However for the parameter values related to the subject 2, the spectrum is 
increasing compared to the linear spectrum given by the parameter values related 
to subject 1. 

Using the new scales — VgT^ and = r„, with velocity 


Do sin(g) 


(15) 


we hnally obtain the complex Ginzburg-Landau equation 



( 16 ) 
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where 


P = 
Qi = 
Q 2 = 


oj'^Dq cos{q) — Dq sin^(g) 
4J3 ^ 


1 r 2 Ai £)2 sin^(|) + — 3A^ + 8 A^ sin^ (!''■' 


u 


ai 


+ 


12 r 2 n + 64iAo sin' 


4/9-1 


Aiz/i ^ SAi^i - 4 ^ 21^1 sin^(I) 


( 17 ) 


2ng 


12 r 2 n + 64iAo sin 


4/gi 


^0 , n • 

7 = j + ^isin (-). 

Dnring the last three decades, the CGL eqnation and its modihed versions have 
drawn tremendous attention. These equations describe a variety of physical phe¬ 
nomena in optical waveguides and hbers, plasmas, Boose Einstein condensation, 
phase transitions, open flow motions, bimolecular dynamics, spatially extended 


non equilibrium systems, etc 


2l| . In the present research work, the CGL equation 


is an equation describing the evolution of modulated hormonal waves in a diffu¬ 
sive coupled pancreatic islet /3-cells model. This result suggests that, the insulin 
propagates within the islet /3-cells using both time and space domains in order to 
regulate glucose level. In another regard, oscillations of insulin secretion which are 
likely caused by intrinsic /3-cell mechanisms generate a spatiotemporal dynamics 
of insulin between cells as modihed by exogenous signals such as hormonal and 
neuronal input. 


We have represented in Fig. 2. the variations of constants P, Qi, Q 2 , 7 and of the 
product PQi with respect to the wave vector q. The parameter values are those 
in the Table 1. It is observed for both subjects that the coefficients are positive 
for small values of the wave vector q. The nonlinearity coefficients have very small 
values. It is observed also that except the dissipative coefficient 7 , all the other co¬ 
efficients decrease with the increasing of the wave vector. It is well known that the 
complex Ginzburg-Landau equation has as the modulational instability criterion 
for the plane waves PiQi -|- P 2 Q 2 > 0 (Pi and P 2 are the real and the imaginary 
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parts of the dispersion coefficient, respectively) for which the plane waves in the 
system are unstable. This relation is known as the Lange and Newell’s criterion 


22| . However, in this work the imaginary part of the dispersion coefficient is equal 


to zero, then the Lange and Newel’ 
Benjamin-Feir instability criterion 


’s criterion reduces to PiQi >0 known as the 


28| . According to this instability criterion, for 


PQi > 0 wave planes in the system are unstable while for PQi < 0, they are 
stable. Since the manner with which hormonal waves propagate in the system 
does not depend of the stability criterion, one can expect to hnd in the system 
spatiotemporal modulated wave solutions for any wave carrier whose wave vector 
is in the positive range of PQi- 


3.2. Nonlinear solution of the equation of motion 

The solutions of nonlinear partial differential equations constitute a crucial factor 
in the progress of nonlinear dynamics and are a key access for the understanding 
of various biological phenomena. Many analytical investigations have been carried 
out to hnd the envelope soliton solutions of these equations which are localized 
waves with particle like behavior i.e., preserving their forms in space or in time or 
both in space and time resulting in spatial, temporal or spatiotemporal solitons, 
respectively 2^ . We assume that the form of the envelope soliton solution of Eq. 
(16) has the form of the one proposed by Pereira and Stenho |^, and Nozaki and 


Bekki 
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i"i(6,r2) = 


Ae‘^ 


(18) 


1 -g ■ 

The real part Fir and the imaginary part Fu of Fi{ii,T 2 ) are given respectively 
by 


Flriil-, a ) — ^ 


e '^ + cos(2/i0)e'^ 
-2(cosh(20) + cos(2/i0))- 


and Fii{ii,T 2 ) =-A 


sin(2/i0)e'^ 


-2(cosh(2(;/)) + cos>{2^(f))) 
(19) 
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where (f) = q^i — 0 JT 2 , ^ = —jS ± ■\/2 + and (5 = —|^- 

Using the solution of Fi given by Eq. (19) and from Eq. (4), we have 

y = 2e{Fir cos 9 — Fu sin 9) + [Fq + 2{F2r cos 29 — F 2 i sin 29)] + O(e^), (20) 

where F2r and F 2 i are the real and imaginary parts of F 2 , respectively such 
that 


F 2 r — Ci{FI^ — F^^) + 2 c 2 FirFii aud F 2 i — C 2 {F‘^- — F'^^) + 2ciFirFii, (21) 


with 

_ Ai - 4D2sin^(|) _ uiui 

“ 3fi§ + WDosin'd) ““ ~ 3fig + 16Dosm‘‘(f)' 


( 22 ) 


Inserting Eqs. (19) and (21) into Eq. (20), we obtain for the insulin dynamics the 
following solution 


yn{t) =eA 


COs(9n — 2Q(0„)e'^" + cos 0^6 


(cosh 2(f)n + cos 2a(f)n) 

+ (ci cos 29n + C 2 sin 29^) x 


Ai 


+ eA^ 


(ci sin 29n — C 2 cos 29n) x 


_1_ 

no(cosh 2(f)n + cos 2a(j). 
2 cos 2a(f)n + cos 

(cosh 2(/)„ + cos 2 q;0 , 

2 sin 2a(pn + sin 
(cosh 2(j)n + cos 2a(j)n 


n) ■ 

'n \ -1 


n 

n \ -1 


(23) 


where = eq{n — Vgt) — ueH. 


In Fig. 3, we have represented the evolution of the solution at different times 
according to the parameter values related to the two hrst healthy subjects of 
Ref. j^. As we observe in this figure, the solution is well a localized modulated 
solution, and its propagates structurally stable. As interestingly remarked in the 
present work, the modulated solution involving in the system appears in the form 
of a breather-like coherent structure and it propagates with the same dynamics for 
the different parameter values related to the two healthy subjects. This assumption 
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can lead to the conclusion that the insulin propagates in pancreatic islet /5-cells 
using localized modulated solitonic waves. 


Let us recall that the localized modulated oscillations obtained in this work are 
involved in many other biophysical systems. Under certain conditions, they can 


move and 
strated in 


ransport energy along the system 


25 


m, 


271]. As recently demon- 


25(1, breathing modes are also responsible of energy sharing between 


a-polypeptide coupled chains. Also, localized oscillations can be precursors of the 
bubbles that appear in thermal denaturation of DNA and they have been shown 
to describe the breaking of the hydrogen bonds linking two bases 2^. It has been 
also shown that these localized oscillations can move along microtubule systems 
[^ . Then, breathers should be understood as triggering signal for the motor 
proteins to start moving as interestingly hnd also in this work. 


In Fig. 4, we have increased the value of e from e = 0.077 to e = 0.099, with the 
same parameter values used in Fig. 3. It clearly reveals in this hgure the influence 
of small perturbation in the dynamics of the hormonal wave. One can easily see 
that for both subjects the solutions still remain the breather excitations. However 
these breathers are now represented by modulated solitons so that the envelopes 
cover less oscillations of the carrier wave as those observed in Fig. 3. It is also 
observed that the amplitude of the wave has increased. 


4. Conclusion 

We have studied the intercellular insulin dynamics in an array of diffusively coupled 
pancreatic islet /5-cells. The model was formulated by a system of discrete ordi¬ 
nary differential equations where the cells are connected via gap junction coupling. 
Motivated with often non-observance clinical effects due to some pathological dis¬ 
eases j^, the work was devoted to derive a clear analytical solution describing 
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the insulin dynamics in pancreatic islet /9-cells. Applying a powerful perturbation 
technique, we have found that the complex Ginzburg-Landau equation is the equa¬ 
tion which describes the insulin dynamics. It has been revealed that the solution of 
the hormonal wave is well a localized modulated solitonic wave called breather. In 
another regard, the breather has been revealed as mechanically important in other 


biophysical systems such as collagen 


251, DNA 261, microtubule 271 systems. The 


correlation with the present work may indicate an important role of breathers and 
other nonlinear excitations in the dynamics of pancreatic islet /9-cells. In a forth¬ 
coming work, we intend to investigate long-range effect, since intercellular waves 
travel also with non contacting cells j3, 301 indicating the long-range interaction 
in the system. 
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Figure 1: (Color online) The dispersion relation of the hormonal wave. Parameter values are: 
D = 0.12 and fig = ai& 2 - Subject 1: ai — 0.0226, 62 = 0.0437. Subject 2: oi = 0.0509, 
62 = 0.2062. 
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Figure 2: (Color online) Variations of coefficients (a) P, (b) Qi, (c) Q 2 , (d) the product PQi 
as a function of the wave vector q of the carrier wave. The parameter values are given in Table 
1. Subject 1: oi = 0.0226, 02 = 3.8 x 10“®, bi = 0.0022, &2 = 0.0437. Subject 2: oi = 0.0509, 
02 = 1.29 X 10"'^, h = 0.0096, 62 = 0.2062. With D = 0.12. 
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(b) 


Figure 3: (Color online) The solution as a function of the position at different times. The 
parameter values are: Subject 1: D = 0.12, q = 0.035, e = 0.077, oi = 0.0226, bi = 0.0022, 02 = 
3.8 X 10-®, 62 = 0.0437. Subject 2: D = 0.12, q = 0.2, e = 0.077, oi = 0.0509, 02 = 1.29 x 10"^ 
bi = 0.0096, 62 = 0.2062. 
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Figure 4: (Color online) Effects of small perturbation on the hormonal modulated wave, e = 0.099 
and t = 5. The parameter values are the same as in Fig.3. 
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